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Abstract 

A variational framework is introduced to describe how a surface bends when it is 
subject to local constraints on its geometry. This framework is applied to describe the 
patterns of a folded sheet of paper. The unstretchability of paper implies a constraint 
on the surface metric; bending is penalized by an energy quadratic in mean curvature. 
The local Lagrange multipliers enforcing the constraint are identified with a conserved 
tangential stress that couples to the extrinsic curvature of the sheet. The framework is 
illustrated by examining the deformation of a flat sheet into a generalized cone. 

1 Introduction 

There is a striking similarity in the patterns we observe in a crumpled ball of paper, a crushed 
aluminium can or, for that matter, the creases that develop on the back of one's shirt when 
one's back has been up against the wall a little too long. This is not an accident. 

The relevant property that these materials share is that they offer resistance to tangential strain. 
This is captured geometrically in the statement that the metric, which determines distances 
between points on the surface (and thus is a quantity intrinsic to it), does not change under 
deformation: the only permissible deformations are isometries [T] [31 II] ■ Not only must the 
surface area remain fixed, right angles must remain right on the surface. 

In general, the number of ways that a surface can bend without stretching is limited. In 
particular, if a flat sheet is bent in more than one direction it must stretch somewhere. This is 
a direct consequence of the geometrical fact-Gauss's inspired insight-that the product of the 
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curvatures at any point on the surface depends only on the metric there [5]: two quantities that 
are defined extrinsically are constrained by the intrinsic geometry. In particular, if the surface is 
flat to start with, so that one or both of its curvatures vanishes, it must remain flat. However, 
the number of ways a surface can be flat is itself limited. It is either part of a cylinder, a cone 
or, more generally, a tangent developable surface (think Frank Gehry [6]), ruled by straight lines 
[UEE]. Any deformation of the surface is constrained locally by the fact that it must coincide 
with one of them. 

The physically deformed surface is, of course, rarely any single one of these elementary flat 
surfaces. This is because the forces applied to the surface usually oblige it to fold along more 
than one direction somewhere. It turns out, however, that the most energy effective way to 
adjust to such forces is by confining the regions where stretching occurs within a series of 
sharp peaks and ridges [S E]. A tangent plane cannot be drawn at these singular points; the 
curvatures will generally diverge. The complete deformed surface will form a quilt of irregularly 
shaped flat patches meeting along the boundaries that these points provide. As described by 
Cerda and Mahadevan they provide the hinges about which the flat patches may flap [10J. Once 
such singular points are established, however, irreversible damage will be done to the surface: 
iron out the creases we may but the paper will remain scored, the can scarred and, perhaps this 
is stretching the analogy, shirts do become threadbare. 

A theory which describes the detailed internal structure of this network of peaks and ridges will 
be complicated and it will depend on the detailed material properties of paper. An approxi- 
mation, reliable for small strains, is provided by the Foppl-von Karman equations The 
length scale of these regions is on the order of the thickness of the sheet; viewed from outside, 
however, this internal structure is irrelevant: peaks become points and ridges can be treated as 
curves along which boundary conditions are set on the flat surface; the only surviving material 
parameter in the physical description of the sheet is a modulus characterizing its stiffness. The 
important point is that on these longer length scales the only relevant degrees of freedom are 
the geometrical ones; how the surface bends becomes an essentially geometrical problem, a fact 
that tends to get lost when it is treated from a continuum mechanical point of view. In this 
paper we would like to explore this geometrical approach to the problem. 

The bending energy of a surface is deceptively simple; modulo a topological energy there is 
a single quadratic in curvature associated with a surface [121 [T3]. However, determining the 
bending patterns of a surface when bending alone is penalized provides a subtle non-linear 
problem in itself and there are few useful analytical results. The problem we are interested in 
addressing has an added complication: there are local geometrical constraints on the surface; 
the metric is fixed. 

When the constraints are global it is straightforward to enforce them in the minimization of the 
energy by introducing Lagrange multipliers. If the enclosed volume and the area are prescribed, 
these multipliers are identified with the equilibrium pressure and tension: they are global. This 
is a situation that is familiar (in the context of fluid membranes see, for example, [H]). With 
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local constraints these multipliers get replaced by local fields with their own conservation laws 
reflecting the constraint. 

In this paper we will set up a variational framework to describe the equilibrium shapes assumed 
by a surface subject to local constraints. In equilibrium, we will show that the multipliers are 
identified with a conserved tangential surface stress. The constraints fixing the metric thus set 
up tangential stresses within the sheet. These add to the stresses associated with bending. 
They couple linearly to the extrinsic curvature in the shape equation describing the equilibrium 
configuration. 

We will illustrate the framework by examining the deformation of a flat sheet into a generalized 
cone [9]. This is the simplest flat deformation of a plane exhibiting a localized distribution 
of energy. Our focus will be on the stresses set up in the cone rather than the geometrical 
details of the configuration itself. We will examine in detail how the different contributions to 
the stress conspire to transmit force and torque along closed curves on the cone. In particular, 
we will show how our framework can be applied to two physically interesting setups that have 
been discussed in the literature. The first of these considers the depression of a circular sheet 
into a circular frame by application of a point force to its center [TU]. Neither the sheet nor 
the frame needs to be circular. We will show that the consistent coupling to curvature is 
enough to completely determine the radial behavior of the stress tangent to a sphere centered 
on the apex of the cone. The radial behavior of the remaining projections of the stress tensor is 
then completely determined by the conservation law. We will show that additional off-diagonal 
stresses, that have not been studied previously, are consistent with the conical geometry. The 
second example considers the 'draping' patterns assumed by a sheet supported at a point that 
falls under the influence of gravity [15]. Additional stress fields are now needed to counteract 
the effects of gravity which would tend to unflatten the disc. We show, however, that they 
are also completely determined by the consistent coupling to curvature. Remarkably, the shape 
equation we obtain assumes a universal form, with or without gravity; the physical details of 
the problem enter the equation only through the stress tangent to the sphere associated with 
the constraint. In particular, with gravity, this equation differs from its (technically inaccurate) 
counterpart obtained by pre-averaging over the radial direction. 

2 Bending a constrained sheet 

We are interested in determining the configurations that minimize the bending energy of an ini- 
tially flat unstretchable surface under the influence of some set of external forces or constraints. 
This might be gravity, forces that act locally or the confinement of the surface within a fixed 
volume. We will first examine the problem without the complications introduced by gravity. 
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The bending energy associated with the configuration X is given by 

H[K} = 1 -JdAK\ (1) 

where K = C\ + C 2 is the sum of the principal curvatures (twice the mean curvature), and dA 
is the area element induced on the surface (our notation is summarized in [IS]). We have set 
the rigidity modulus to be one. 

Bending must be an isometry: not only does the geometry of the sheet remain flat almost 
everywhere when it is bent, it also resists shear and stretching. Thus it is not sufficient to 
demand that the Gaussian curvature vanishes; one needs to fix the metric itself. The bending 
energy H[X] defined by Eq. (JT]) needs to be minimized subject to this constraint. The frame- 
work will also need to accommodate discontinuities and singularities which arise due to global 
obstructions on this requirement. 

Distances on the surface are described by the metric, g ab = e a ■ e&, where e a = <9 a X, a = 1,2, 
are the two tangent vectors adopted to the particular parametrization of the surface. In a sheet 
of paper the metric is fixed, g^ b say. We thus construct a functional to reflect this constraint 

H C [X] = H[X] -\J dAT ab (u a ) (g ab - g$) . (2) 

Here T ab {u a ) is a set of local Lagrange multipliers implementing the constraint on the metric. 
Both g ab and g^ as well as the multipliers T ab form tensors, so that the trace appearing in 
Eq. (P3 is a scalar. Thus, by construction, the functional H c is independent of the parametriza- 
tion. 

With the constraint in place, it is possible to introduce a deformation X — > X + 5X in the 
shape without needing to worry if the metric is behaving. 

In the absence of external sources of stress, the constrained equilibrium may be expressed as 
the conservation law [IT] 

V a P = 0, (3) 
where the stress tensor f a is a sum of two terms: 

f a = f a + T ab e b . (4) 

The bending contribution is given by (for our notation see [16J) 

f a = K(K ab - l -g ab K) e b - {V a K) n . (5) 

There is a particularly direct way to derive these equations. It involves a natural extension of 
an approach to the equilibrium of fluid membranes using a set of auxiliary variables introduced 
by one of the authors [T8] , 
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The constraint adds a tangential stress proportional to the multiplier T ; tension will be 
introduced in the surface whenever its action is antagonistic to bending. This tension will 
generally be inhomogeneous and anisotropic. 

The 'shape equation' is given by the normal projection of the conservation law (T3]): 

S = S -K ab T ab = 0, (6) 

where 

S = - V 2 K + \K{K 2 - 2K ab K ab ) . (7) 

Here V 2 is the Laplacian on the surface constructed with the induced metric. The stress T ab 
couples linearly to the extrinsic curvature in Eq. (|6]). 

The tangential component of the conservation law (T3D is the statement that T ab itself is con- 
served, 

V a T ab = . (8) 

The bending energy does not enter this equation. It is an intrinsic statement involving the 
intrinsic geometry through the covariant derivative. We note that T ab involves three degrees 
of freedom; this is the correct number of components to pin-down the three components of 
the metric. From a field theoretical point of view it is also possible to view the two con- 
straints provided by Eq. ([H]) abstractly as a requirement for consistency on the T ab implied by 
reparametrization invariance [17J. Note that the tangential part of f defined by Eq. ([5]), unlike 
T ab , is not conserved. 

It is worth pointing out that this variational approach to bending with a local constraint is more 
direct than its continuum mechanical counterpart. One does not need to invoke any adapted 
system of coordinates. One also sidesteps the need to introduce constitutive relations; once the 
energy function is given, the relationship between stress and strain follows. 

One now needs to solve the coupled system of equations §5§ and ([H]) subject to appropriate 
boundary conditions. Eq. (E]) will pick up a source if external forces are operating on the sheet. 
In addition, the geometry will typically display geometrical singularities where Eq. (|8]) will pick 
up a distributional source. We will apply this framework to the conical geometry in the next 
section. 

It is also instructive to examine how this setup would be different had we fixed only the Gaussian 
curvature. This is addressed in appendix[A] The relevant extension of our framework to handle 
the situation in which the initial surface is not flat is a straightforward generalization. 
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3 Folding a cone 



3.1 Geometry of a Developable Cone 

We begin with a brief review of the conical geometry described by Ben Amar and Pomeau [5]. 
For further details see [TP] . 

A cone possesses a natural parametric representation in terms of a closed curve r : s — > u(s) 
on the unit sphere. Let s be arc-length along this curve and r be the distance from the origin 
along the ray pointing in the direction u(s) (see Fig. [lj. The image of the mapping 

(r, s) — > X(r, s) = ru(s) (9) 

describes a cone. We will suppose that this cone is the result of an isometric deformation 
of a planar disc (not necessarily circular) of radius R(s). The length of T is then fixed at 
L = 2n. This means that the closed curve which represents the cone lives always within a 
single hemisphere. A flat planar disc is represented by a great circle. These are the only circular 
configurations consistent with the constraint on the length. If the disc is circular of radius R 
and the apex coincides with its center, the maximum value of r is given by R. 

The tangent vectors to the cone adapted to this parametrization are u itself and t = u', where 
prime denotes a derivative with respect to s. Note that u ■ u' = so that the induced metric 
on the surface is specified by the flat line element 

dl 2 = dr 2 + r 2 ds 2 . (10) 

Let n = u x t denote the normal to the surface. The extrinsic curvature tensor is given by 

K- = r(J»), (11) 

where 

K = -n • t' = -(u x t) • t' = -u ■ (t x t') . (12) 

The flat direction is along u. The Gaussian curvature K G = det (K%) vanishes. The curvature 
K = n/r diverges at the apex of the cone. This translates into a logarithmic divergence in the 
bending energy at this point. 

The three vectors u, t and n form a right-handed basis. Given k, the curve can be reconstructed 
using the structure equations 

u' = t (13a) 
t' = -Kn-u (13b) 
n = «t. (13c) 
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Note that the vectors u and n are in general not equal to the normal and binormal vector of the 
curve [7]. Thus k should not be mistaken for the Frenet curvature, Kp. It is simple, however, 
to show that the two are related by up = k 2 + 1. 



3.2 The cone as a constrained minimum 



We have seen that cones can be identified with curves on a unit sphere. As described in 
appendix [Bl by integrating over the radial direction H becomes a functional of curves on this 
sphere. As such, one could vary the Hamiltonian with respect to these curves to obtain an 
Euler-Lagrange equation which is satisfied by the curves which minimize the energy. This curve 
is then identified with the equilibrium shape of the conical surface. While this approach yields 
the correct shape equation in this case, it does not always provide the correct answer. As we 
will see, variation of the Hamiltonian and restriction of the geometry (integrating out the radial 
dependence in the conical geometry) do not commute in general even when the disc is circular. 
Even in this case, however, it does not possess the scope to provide a complete analysis of the 
problem. 

We will instead apply our general framework to study the conical equilibrium. We first confirm 
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that cones are in fact equilibrium geometries. It is straightforward to show that 

S = - V 2 K + l -K(K 2 - 2K ab K ab ) = -1 L" + l - k 3 + . (14) 

Note the appearance of the term linear in k originating in the radial part of the Laplacian. The 
shape equation thus reads 

k" + -k 3 + k) - nrT ss = 0. (15) 



r 3 V 2 

Only the stress tangential to the generating curve on the unit sphere contributes 



The shape equation describing the conical equilibrium must hold for all r. Indeed, consistency 



requires that r 4 T ss depends only on s. Let us write r 4 T ss = — Ci(s) 



At this point it is also convenient to replace the components of the stress tensor T ab by the 
corresponding scalar quantities representing its projections along the directions tangential to 
the sphere and along its radius, t and u: T\\ = t a t b T ab , % = t a u b T ab , and T x = u a u b T ab E 
In particular, the tangential stress 

r„ = t s t s T ss = r 2 T ss = -C|| (a) 1 (16) 

diverges as r~ 2 . We see immediately that the consistency of the shape equation with a conical 
geometry places a very strong constraint on the admissible form of Ty, determining completely 
its radial dependence. 

Let us now examine the conservation law for T ab to see what it tells us about Tli, T± and Tim. 
The tangential and radial projections of Eq. flS]) read: 



1 ™, d „ 2 
- T ll + ^ + - 

dT, 1 ,„ „ , 1 



-T,; + — T||± + -T|| ± = 0; (17) 



^r--( T ^- T ii) + - T iU = o- ( 18 ) 

Often the symmetry of the problem will imply Tim = so that, modulo Eq. (TT7I) . Cn is a 
constant independent of s. If the circular frame discussed in [10] is replaced by one that is 
irregular, however, the balance of forces on the frame will require Tin ^ 0. 

Remarkably, when C\\ is constant, as it will if the disc is circular, the shape equation assumes 
exactly the same form as the Euler-Lagrange equation for planar Euler Elastica with tension 
cr oc Cm — 1 and k in place of the Frenet curvature, k f . 



1 With respect to the adapted coordinate system, t a — (0,r x ) and u a — (1,0) so that t a — (0,r) and 

U a = (1,0). 
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The Euler-Lagrange equation for planar Euler Elastica is completely integrable for k allowing the 
cone to be reconstructed using the structure equations. The tension in this fictitious Elastica 
can be either positive or negative, depending on the boundary conditions. This should be 
compared to the results of Cerda and Mahadevan who do not admit negative values [10J. Note 
that if k f is used in place of k the equation is different. The shape does not correspond to 
Euler Elastica on a sphere 



More generally, the off-diagonal constraining stress is given by 



. ylnr C l[± (s) 



(19) 



where C\\± is some other function of s. In a moment, we will see that if Cu±(s) is non-vanishing, 
it will play a role in the balance of torques on the cone. 



The second conservation law ( JT8f ) identifies the radial scalar T± as a sum of five terms 

'lnr 1\ C||(s) C„ 



CAs)" 



In particular, if Cm is constant and C\\ 



0, 



C|| , Cj.(s) 



(20) 



(21) 



The r~ 2 behavior term is slaved to Tji . This part of the stress tensor T ab is traceless: g ab T a b = 0, 
a property that can be tracked back to the scale invariance of the bending energy. The r" 1 
behavior indicates the existence of a source term associated with external forces. 



It is worth emphasizing the following points: the stress associated with the constraint on the 
metric is explicitly anisotropic and inhomogeneous. It is not necessarily diagonal with respect 
to the orthonormal frame adapted to the cone. While it is not monotonic, it does vanish 
asymptotically. 

It is also useful to possess explicit expressions for the surface components of the stress, fy = t a f a 
and fj_ = u a f a representing the force per unit length transmitted across curves of constant s 
and r respectively |2D]- A straightforward calculation gives 



Q k 2 + r 2 T^j t + r 2 T\\ L u-k'u 



■- K 

2 



r 2 T ± 



u + r 2 7|, ± t 



k n 



(22) 
(23) 



As a check of consistency, note that the divergence may be expressed as 



v Q p = -f,; 



r 



Of i 1 

~Q ^ _tj 

or r 



(24) 



In particular, it is simple to confirm that the projections of Eq. (1241) along n, t and u reproduce 
Eqs. (US]), (HZ]) and (HHD respectively. 



4 Force and torque balance at constant r 

Now look at the balance of forces along a closed curve T r of constant r. The physical arc-length 
is dr = rds. The normal to this contour pointing away from the apex is the vector u. Then 
the total force transmitted across the curve is given byo 

F = <j) dr fj_ = — (j) ds — -k 2 u*j + r <j> ds T± u + r <j) ds T\\± t 

= — (b ds dsT ± u + r j> ds T\\± t , (25) 

where we have used Eq. (123]) on the first line and Eq. ( |13 b[) on the second. A much more 
transparent expression is possible in terms of the solution of the conservation law Eq. ([8]) for 
T ab . Look at the integral along any closed curve of the Euler-Lagrange derivative £, weighted 
by n. In equilibrium, 

ds ( k" + - k 3 + k + nr 2 T\\ J n = . (26) 



2 

Now, performing an integration by parts on the first term and using Eq. (I13cj) . we obtain 



ds k" n = - & ds KK't . (27) 



Furthermore, using Eq. ( ]13bh to express n as a linear combination of u and t' , we obtain 

ds Q k 3 + k + Kr 2 T^j n = - j> ds Q k 2 + 1 + r 2 T^j (f + u) 

ds KK't + r 2 <j) ds Tjj t 

- j ds Q k 2 + 1 + r 2 T\^j u . (28) 
Thus, the tangential stress is constrained as follows on a closed cone: 

ds U + T^ 2 ) u + r2 ldsT\\u-r 2 *dsT|jt = 0. (29) 



2 Any closed curve homotopic to r r can be used to determine the force replacing i± by l a f a , where 1 = l a e a 
is the normal to the curve on the cone pointing towards increasing r. 
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Modulo this identity we can finally write the total force in terms of the function C± appearing 
in Eq. {23]) as 

F = r I ds [(T ± + T,|)u + (% - Tjj) t] = j ds C ± (s) u , (30) 

where we have used the expressions for the projections of the stress tensor ( JTBT ). ( JT9l ) and ( I20l) 
consistent with the shape equation and the conservation of T ab . The line integration does not 
depend on r. If an external force is acting, the term proportional to r^ 1 in T± does not vanish. 
If a force F is applied at the apex of a cone supported on a frame, C±(s) is the distribution of 
forces on the rim of this frame needed to counterbalance it. Notice that the tangential and the 
radial projections of the conservation law can be written as 

V a T ab = [C ± (s) u b + C ll± (s)t b ] ^ . (31) 

In exactly the same way, we can consider the torque. The conserved torque tensor is given by 
[IZII22| 

m a = X x f a + Kg ab e b x n (32) 
and the total torque about the origin can be written as: 



M = f dr u a m a . (33) 



Thus 

M 



ds (ru x fj_) + Ku x n 
-2 (b dsnt + r 2 <j) dsT\\ ± n 



I dsC ll± (s)n. (34) 



The integrated contribution from fg vanishes on account of the structure Eq. ( I13cp and the 



closure of T r . Furthermore, the term proportional to lnr/r 2 in T\\± [see Eq. ( IT^l) ] does not 
contribute to the torque as one can show by integrating the Euler-Lagrange derivative (TTBT) 
weighted by t along a closed curve. 

Eq. (!34|) implies that the cone geometry may provide an equilibrium when an external torque 
is applied. In that case, the off-diagonal stresses do not vanish and the value of the torque is 
captured by C\\±(s). 
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5 Brought down by gravity 

Suppose that the sheet is held at a point and allowed to fall under gravity. The bending energy 
is minimized by the planar disc; gravitational potential energy is minimized by a sheet that 
hangs vertically; as described in [15], the competition between the two gives rise to draping 
patterns. 

If the weight of the sheet is taken into account, its gravitational potential energy is given by 

H g [X] = -pg J dAX-z, (35) 

where p is the density of the sheet and g is the acceleration due of gravity. To counterbalance 
the weight of the sheet an external force F ext in the z-direction has to be applied (see Fig. [T]). 
This force acts at the apex of the cone (see previous section). 

Under the deformation, X — > X + 5X, the variation of H g is given by 

SH g = -pg J dA (KX ■ z + n • z) n • SX . (36) 

The effect is to add a normal source to the conservation law ([3]): 

V„f a = pg(KX- z + n z)n 

= pg(nu ■ z + n • z) n = j(s) n . (37) 

The Euler-Lagrange equation (TT5!) is replaced by the equation 

-^(- /, + ^ 3 + -)-7^ll=^)- (38) 
This differs from the result that Cerda et al. have derived. They average over the radial direction 

na. 

As before, consistency will place a constraint on Tji. However, we now require thalEI 

SiCf) _i r = Tr o_ i r (39) 

T K ft 

Even in the regime where the analysis provided in appendix [B] is valid, with gravity we begin 
to see how such an approach breaks down. The pre-averaging over r provides a different (and 
incorrect) shape equation. 



5 ln the following, all variables with an upper index g = refer to the case without gravity. 
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The conservation law fl£ 
However, because of Eq. 



for T ab is unchanged as gravity only contributes a normal source, 
in place of Eqs. ( Tl§l) and (I2"01 we have 



r 

3' 



and 



The force transmitted across a curve of constant r is then given by 






~] 


1 


(k, 


II- 

) 


r 




K 


^3 






2 



(40) 



(41) 



■ tp ds 



JT3bl 



dr fj_ 

ds Ci (s) u + r 2 



(T ± + ^f =0 )u+(7|| ± -^ g=0 )t 

ds 



/C 3 \K 

ds C ±(s) u + — (D ds j n — — ® ds 



u 



' 1 ' 
] - \ u 

K 



dAjn,(42) 



where E r is that part of the surface which is enclosed by the curve T r . Note that the line 
integration does now depend on r but for good reason: the surface below the curve T r exerts 
the force F on T r due to gravity. Its absolute value decreases quadratically for increasing r and 
vanishes for r = R as expected. 



The torque about the origin is now given by 



M = r 2 i dsT ll± n 



M g=o , r 2 



ds ' 



n = M g= 



dA rjt 



(43) 



The last term due to gravity can also be written as J E gL4 [X x (jn)], again exactly as one 
would expect. 



6 Discussion 

In this paper we have introduced a geometrical framework to examine the bending of an (un- 
stretchable) sheet of paper. We have illustrated the viability of this framework by examining 
the conical deformations of a planar sheet. In particular, we have examined in some detail the 
distribution of stresses within the sheet associated both with bending and with the constraint 
and how they conspire to transmit the external forces and torques acting on the sheet. 

Various characteristic patterns are associated with this geometry: the radial dependence of the 
stress is constrained strongly by the geometry; how it varies along the spherical generating curves 
will depend on the specific boundary conditions associated with the external forces acting on the 
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sheet. We have also shown how the presence of bulk forces such as gravity are accommodated 
within this framework. In the latter case, we point out the pitfalls of prematurely constraining 
the geometry within the variational principle. 

We have not attempted to discuss the physical details within the ridges and peaks where the 
elastic approximation breaks down. Various elements of our framework are, however, likely to 
play a role in any refinement of the model which does address this physics. 

We have limited our discussion of the application of our framework to conical shapes; this is 
because they are relatively simple not because they are the only interesting configurations. It 
would be interesting (and it should be straightforward) to examine the Mobius strip geometry 
described recently by Starostin and van der Heijden within our framework [23] , 

It is appropriate to think of a cone as a kind of elementary deformation of a planar sheet in 
which the energy is localized. A generic patch of flat surface, however, will be described by 
a tangent developable surface. Such a surface is singular along a certain curve: its edge of 
regression which will itself generally exhibit singularities [SIT]. The flat directions are generated 
by the tangent to this curve. Typically, we do not see this curve when we fold paper. The patch 
gets truncated by a ridge or we run off the sheet before the corresponding generating curve is 
reached. However, the folded sheet is, in principal, completely described by this set of curves. 
This description will be the subject of a subsequent publication |24|. 
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A Fixing Gaussian curvature 

Let us examine the consequences of fixing the Gaussian curvature instead of the metric. We 
thus construct the constrained functional 



where Kq is the Gaussian curvature [T(5]- Here W(u a ) is a local Lagrange multiplier constraining 
Kg to vanish. Analogous functional can be defined for a constraint on any other geometrical 
scalar. 





(44) 
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The constrained equilibrium may again be expressed as a conservation law of the form ([3]). The 
stress tensor f a is again given by an expression of the form (@J with an additional tangential 
stress T ab . However, this stress is now determined completely by the multiplier W, 

T ab = _ (yaybyy _ g ab V 2 w ^ _ ^ 

The constraint again introduces tension in the surface but, this time, with a single degree of 
freedom. 

It is appropriate to point out that there is an ambiguity inherent in the definition of the stress 
tensor f a . This ambiguity is due, in part, to the Gauss-Codazzi integrability condition which 
connects the intrinsic and extrinsic geometry [21]. If the variations are treated intrinsically, 
one obtains the expression we have written down for the stress; if, however, they are treated 
extrinsically an apparently very different answer is obtained. The two are entirely consistent; 
they differ only by a conserved null stress that does not transmit forces. The relevant calculations 
have appeared in a different context (see, for example, [18J ) . 

Whereas the conservation of T ab is an important constraint when we fix the metric, here it 
is simple to confirm that T ab is automatically conserved if the geometry is flat. This is a 
consequence of the definition of T ab in Eq. (1451) in terms of a potential. In analogy with 
elasticity theory, it is appropriate to think of W as an Airy potential for T ab 

Note that the constraints break the invariance of the two dimensional bending energy with 
respect to conformal transformations of three-dimensional space. In general, the Gaussian 
curvature is not preserved. It is simple to show that under inversion in the origin X — > X/|X| 2 , 
one has [7] 

K G -> \X\ 4 K G - 2(X- n)|X| 2 K+ 4 (X ■ n) 2 . (46) 

Remarkably, the only geometries which remain flat are the cones we considered with apex at the 
origin, and which thus satisfy X ■ n = 0. Given one cone, its inversion in the origin is another 
|21j . Note, however, that curves on the unit sphere are invariant under inversion in the origin. 
Thus, while the physics might be different, the cone remains the same. 

B Cones as trajectories on spheres 

B.l Derivation of the Euler-Lagrange equations 

If we introduce the cutoff r at the apex and integrate over the radial direction, the bending 
energy of a cone formed by a circular disc with its apex at the centre is given by 

H[u] = j dsn 2 , (47) 
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where a is given by 

a = \n(R/r ) . (48) 

The dimensional dependence is contained completely within this logarithm. This behavior is a 
consequence of the scale invariance of H. H is a functional of curves on the unit sphere. 

We have already pointed out in section [321 that if the boundary conditions are not symmetrical, 
the identification with Euler Elastica breaks down. It is, however, useful to see how the shape 
equation emerges from a variational principle which exploits the mapping from cones into 
trajectories on a sphere. Our task then is to identify the trajectories on a sphere minimizing 
the Hamiltonian H given by Eq. (1471) . 

Although the Euler-Lagrange equation itself is simple, its derivation is quite subtle. As we have 
found elsewhere, it is useful to introduce an appropriate set of auxiliary variables associated 
with local geometrical constraints [18]. Consider the functional 



F[u, t, A, A, f] = j ds 



|(u-txt') 2 + ^(u 2 -l) + ^(t 2 -l) + f.(t- U ') 



(49) 



where A(s) and A(s) are two local Lagrange multipliers enforcing the constraints that u is a 
unit vector and that the parameter s is arc-length on the sphere. The area constraint is thus 
implemented by fixing the range of integration on s to the interval [0, 27r]. 

Variation of F with respect to u gives 

f = —{an 2 + A) u + an n , (50) 

The right hand side of Eq. (1501 has been simplified by using the structure equations (TT3I) for the 
curve. Note, in particular, that f is not constant as it would be for a translationally invariant 
Hamiltonian. In the functional F translational invariance is broken both in the energy itself and 
in the constraint on u. 



Variation with respect to t identifies the vector f appearing in Eq. (1501) as 

f = - (A + 2aK 2 ) t - an'n , (51) 

Later, we will provide a physical interpretation of f as an effective force per unit length along 
a ray of fixed s. 

By substituting the expression for f given by Eq. ( T5Tj) into Eq. ( 150"]) we obtain three equations 
describing the equilibrium. The projections along u and along t determine the multipliers A 
and A. The remaining equation along n is the Euler-Lagrange equation. 

Note that 

f-u = 0. (52) 
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Let us first examine the component of the conservation law along t. We note that 

f'-t^O. (53) 
This equation fixes the function A appearing in Eq. (I5T|) up to a constant. We have 

A=c-y K 2 , (54) 

where c is a constant. We thus identify 

f = (^an 2 — c^j t — an' n. (55) 

Comparing the two stress tensor projections ( 1221 and ( 1231 from Sec. 13.21 with the Lagrange 
multiplier f yields the remarkably simple relation f = ar 2 fy if Tin = 0. This confirms that, 
when this approximation is valid, f is the tangent stress along a ray of fixed s. 



To determine A, note that 



Using Eq. ( 1551) . it follows that 



A = -^ 2 -c. (57) 

Note that this variable does not appear in the final shape equation which involves normal 
projections. For the record, we have 

f' = -^ K 2 - c J u + ao. (58) 
Its normal projection together with Eq. (1551) provides the Euler-Lagrange equation 

K," + ^k 3 + k^ +CK = 0. (59) 

The constant c is a tension associated with the fixed area constraint which is implied by fixing the 
arc-length. If this constraint is relaxed, c = 0; however, there persists a 'tension' proportional 
to a associated with constraining the curve trajectory to a sphere. Earlier we identified the 
origin of this term in the r dependence of V 2 K on the cone (see Sec. 13.21) . We reproduce 
the shape equation ( 1591 obtained earlier if we identify Cm with the ratio of the constants (c/a) 
(when Cm is constant). 

Generally, however, it appears to be beyond the scope of the variational analysis based on the 
reduced Hamiltonian to provide a consistent description of the cone. In particular, it does not 
appear to be possible to accommodate Tim ^ 0. Even when Tim = 0, there is no obvious way, 
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within this simplified framework, to analyze the r dependent stresses set up in the cone as they, 
for instance, show up in the presence of gravity (see section [5]). 



We note that there is an integrability condition associated with Eq. (1551) . The identity § dsf 
associated with closure implies that the integrated source of f must also vanish. Thus 



ds 



(k 2 - 2c/ a) u- 2k n =0. (60) 



This identity can be simplified further using Eq. H13 b[) to give 



j ds (k 2 - 2a /a) u = 0. (61) 

There is no obvious physical interpretation of this identity within this one-dimensional frame- 
work. It is instructive to compare Eq. (loTI) with its analogue Eq. (I2U1) within the general 
framework which we used to place a constraint on the stresses set up in the cone. 

B.2 Rotations and Conservation laws 

Fixing the apex of the cone at the origin breaks translational invariance. The residual symmetry 
is three-dimensional rotational invariance. The shape equation describing the sheet may be 
identified as the conservation law associated with this symmetry. In equilibrium, we must have 

/d 
ds— [f • <fu + a«(u x t) • 6t] . (62) 
ds 

In particular, 5F = under rotation. With 5u = h x u and 6t — b x t, we have 

b • [u x f + an u] (63) 

is constant. Thus 

uxf + <mu = J, (64) 
where J is a constant vector which can be identified as the tangential projection of the torque 



tensor 

f K 

t a m a = X x f,| + Ku = X x — + -u = — J , (65) 



ar 2 r ar 



ifT|, ± = 0. 



Using the structure equations (TT31 and the Euler-Lagrange equation ( l59j) one can easily show 
that J is indeed conserved. Squaring, using the identify f • u = 0, we find 

f 2 + a 2 K 2 = J 2 . (66) 

This is the first integral of the Euler-Lagrange equation. It can be rewritten 

1 T2 _ 2 

(k') 2 + \^ - V = , (67) 

and integrated after a separation of variables. 
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